arXiv:hep-ph/0211398v2 30 Jan 2003 


hep-ph/0211398, MSUHEP-21125, SMU-02/05 


Diphoton production in gluon fusion at small 
transverse momentum 


P. M. Nadolsky*^ and C. R. Schmidt'^ 

^Department of Physics, Southern Methodist University, Dallas, TX 75275, U.S.A. 

^Department of Physics & Astronomy, Michigan State University, East Lansing, 

MI 48824, U.S.A. 


Abstract 

We discuss the production of photon pairs in gluon-gluon scattering in the context of 
the position-space resummation formalism at small transverse momentum. We derive 
the remaining unknown coefficients that arise at 0{as), as well as the remaining 
0(ag) coefficient that occurs in the Sudakov factor. We comment on the impact 
of these coefficients on the normalization and shape of the resnmmed transverse 
momentum distribution of photon pairs, which comprise an important background 
to Higgs boson production at the LHC. 


The Standard Model production of photon pairs with a large invariant mass 
plays a vital role in physics studies at the Large Hadron Collider (LHC). It pro¬ 
vides a large background to the production of Higgs bosons, where the Higgs 
boson subsequently decays in the diphoton channel (pp —^ HX 11^) ■ 
Despite the small branching ratio of the Higgs boson to two photons, this 
mode is the most important one for Mh < 140 GeV, due to the narrow width 
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of the Higgs boson and the fine mass resolntion of photon pairs in the LHC 
detectors [1], which allow a Higgs boson peak to be found above the contin¬ 
uum background. The efficient discrimination of Higgs boson events from the 
background relies on the accurate knowledge of the kinematic distributions of 
both signal and background. In a recent paper [2], we and our collaborators 
discussed the diphoton background and calculated the transverse momentum 
distribution of the photon pairs in the framework of the Collins-Soper-Sterman 
(CSS) resummation formalism [3,4]. This resummation is necessary to handle 
correctly the large effects of soft and collinear QCD radiation at diphoton 
transverse momenta Qt of about Mh/2 or less. 

In Ref. [2], significant attention was paid to the production of photon pairs 
in gluon-gluon fusion gg 77 X. This subprocess first arises at C1 (q;|) in 
the perturbative expansion in the QCD coupling. Thus, it is formally of a 
higher order than the quark annihilation subprocess qq 77 X, which enters 
at Oi^oPg). Despite the extra factors of 0:5, the two contributions are compa¬ 
rable numerically, because of the large gluon luminosity in the relevant mass 
range at the LHC. Furthermore, the lowest order (LO) gg 77 contribution 
occurs through a one loop box diagram, which is infrared finite and is not 
related through factorization to the 0{a^g) and 0{ag) diagrams in the quark 
annihilation channel. Therefore, it can be treated as the LO diagram of an 
independent perturbative contribution to diphoton production. 

Recently, the complete next-to-leading (NLO) cross section for the gluon fusion 
subprocess has been calculated [5]. That calculation utilized the cross sections 
for the 0{a‘^ag) real emission subprocess gg ■y'yg [ 2 , 6 ] and the recently- 
computed two-loop virtual corrections to the Cl(a|) box diagram [7]. In this 
paper, we use the results of the above publications to derive all the NLO 
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coefficients in the resnmmed cross section and the remaining nnknown NNLO 
coefficient in the perturbative Sudakov factor. 


In the CSS formalism, the gluon-fusion cross section at small transverse mo¬ 
mentum can be expressed as a Fourier-Bessel transform of a form factor, 
W{b,Q,v,XA,XB), in terms of the impact parameter b: 


d(T{gg 77 ^) 
dQ^dydvdQl 


d% 


(&, Q, V, xa, xb). 


( 1 ) 


The perturbative part of W(b,Q,v,XA,XB) can be written as 


W{b,Q,v,XA,XB) = - ^ - 

r rClQ-dli^\.. , ,,, ClQ^ , 
xexpl- — A{as{y))\n—^ + B{as{y)) \ 

^ ® oisi^p)) [Cg/b ® fb] {xb, b; Up, as{yp)). (2) 

a,b=g,'E 

Here Q, y, and Qp are the invariant mass, rapidity and transverse momentum 
of the photon pair, respectively; S is the square of the pp center-of-mass (c.m.) 
energy; n = (1 — cos6**)/2, where 9* is the polar angle of one of the photons in 
the 77 c.m. frame; and xa,b = jVS- The momentum scale at which the 

QCD coupling as is evaluated is shown explicitly in each of the terms. The 
convolution is dehned in the conventional manner, 

[f®g]ix)= [ ^fiOaix/^) . ( 3 ) 

The summation over the indices a and b goes over the gluon parton distribution 
function (PDF) fg{x,p,p) and the quark singlet PDF f-E{x,p,p), which are 
evaluated at a momentum scale pp. The parameters Ci and C 2 in the Sudakov 
term are constants of order unity. In general the scales po and pp should be 
of order Q and 1/b, respectively, so as not to introduce large logarithms in 
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eq. (2). For the process gg 77 X, the normalization factor is 


O-Q 


647r(5^ 


( 4 ) 


where qi are the charges (in units of e) of the quarks that run in the box loop, 
and the second summation is over the helicities Ai, A 2 , A 3 , A 4 of the gluons and 
photons. The LO helicity amplitudes are given in Eq. (3.15) 

of Ref. [7]. They can be expressed as functions of n = (1 — cos6**)/2 = —t/s, 
where t and s are Mandelstam variables of the LO 2-to-2 process. 


The functions 71 (q;s'), B{as), and Cgjaix, b; fip, as) can be expanded as a per¬ 
turbation series in 0 : 5 : 71 ( 0 : 5 ) = (as/ 7 r)” 7 l(”\ I3{as) = (as/Tr)"" 

and Cg/a{x,b-, fiF,as) = 5ag5(l - x) + En=i («s/7r)”C^/^(a;). For brevity, we 
suppress the explicit dependence of Cg^'l{x) on b and jj^p. The coefficients 
and in the Sudakov factor have been known for some time [9]: 


=-/3o - 2iVOn 


b,C2 

Cl ’ 


^( 2 ) = jv W, - —iV/ - /3o In ^ 

VV 36 127 18 ^ Cl 


(5) 

( 6 ) 


where Nf is the number of active quark flavors, Nc = 3, Cp = 4/3, Pq = 
(lliVc — 2iVj)/6, and 60 = = 1.2292... . We hnd that the 0 {as/'K) 

convolution functions C^/a(a^) can be written as 






= 5{l-x) ("^^k^ + ^^olnl-Z^oln^-Wln^^' 


-P];j(a;)ln 


= NfCpx - Pl)’-p{x) In 


Upb 

fi) 


(7) 

( 8 ) 


where Pg^l{x) and Pg^hix) are the 0{as) splitting functions. All terms on 
the right hand side of Eqs. (7) and ( 8 ), except for Vgg^^^{v), can be obtained 
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from the order-by-order independence of the fnnction W(b, Q, v, xa, xb) on the 
parameters /io, and C 2 , as well as the nniversality of the off-diagonal 

contribution particular, the term /9oln(/io/Q) occurs because the 

LO cross section is 0{ag), and it implies that the natural scale for evaluating 
0:5 in Eq. (4) is /xq = Q- The function Vgg^^^{v) can be obtained from the 
two-loop corrections to the gg 77 matrix element of Ref. [7]. We find 

277eE 


V, 


99^77 


= + 




(9) 


where the summation is over the helicities Ai, A 2 , A 3 , A 4 of the gluons and 
photons. The helicity amplitudes and 

FSL ^ F_fj^ 5 v 2 A 3 A 4 given explicitly in Eqs. (3.15, 4.7-4.16) of Ref. [7]. 


In previous studies [2,8], before the diphoton two-loop virtual corrections were 
available, the functions for the process gg 77 X were approximated 

by their counterparts for Higgs boson production, gg HX, calculated in the 
mtop —> 00 limit. The rationale for this was that both processes are initiated 
by a gg initial state and occur through a quark loop at LO. Thus, the NLO 
corrections were expected to be comparable. The functions for Higgs 

boson production are also given by Eqs. (7) and ( 8 ), except for the replacement 
of Vgg^^^ by [ 10 ] 


Vgg^B — 5iVc — 3Cf XcTt'^ — 11 -f Stt^. 


( 10 ) 


Clearly, the use of the Higgs C-functions would be justified if Vgg^n is numer¬ 
ically close to Vgg^^^. To estimate the validity of this approximation, we plot 
in Eig. 1(a) the quantities Vgg^.yy/A and Vgg^H/^ as functions of the variable 
n = (1 — cos9*)/2. Eor the “canonical” choice of parameters Ci = 60 , <^2 = 1, 
/io = Q, and = bo/b we have C^/^(a;) = <5(1 — a;)V/4; hence the magnitude 
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of V completely determines the size of the gfgf-initiated NLO correction. 

Fig. 1(a) shows that is symmetric with respect to n 1 — n and 

becomes singnlar in the limits n —0 and n —1. These singnlarities, which 
are proportional to powers of Inn, do not contribute to the experimental cross 
section; they are removed by cuts on the transverse momenta of the observed 
photons 7 i and 72 . For instance, the selection cuts used in Ref. [2] were > 
25 GeV. At LO this imposes the constraint (1 — i?)/2<n< (l + i?)/ 2 , with 
R = {1 — • The excluded regions for Q = 120 GeV are shown 

by the shaded areas in Fig. 1(a). We see that in most of the allowed region 
the function is nearly flat, with a numerical value of about 6.65. 

For comparison, we also plot in this hgure which has a value of 

~ 10.15. Thus, the approximation of substituting the Cg/g(a^) coefficient from 
Higgs production overestimates by about 50%. On the other hand, we note 
that the contribution ll/d to V^g^n/d comes entirely from the short-distance 
renormalization to the effective Hgg operator, which has no counterpart in the 
gg 77 process. If we remove this short-distance contribution from Vgg^n/d, 
we are left with W^r^/d 7.d0, which only overestimates by about 10%. 

It is interesting to note that the comparable corrections to the process qq 77 
are considerably smaller than for gg 77. In Fig. 1(b) we plot the analogous 
function Vgg^y,y,/d {i.E., thc coefficient of the <5(1 — x) term in which 

was given in Refs. [8,11]. We see that it is equal to 0.5 — 0.6 in most of the 
kinematical region selected by the LHC cuts, which is much less than the value 
of 6.65 that we found for the pp-initiated process. In this hgure we also plot 
the analogous coefficient V^y/d for the Drell-Yan process, which differs from 
Vgg^-y-y/d by loss thuu 5 — 20% over this kinematic range. 
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Figure 1. Comparison of the functions V/4 (a) for gg 77 (solid line) and Higgs 
boson production gg ^ H (dashed line); (b) for qq 77 (solid line) and the 
Drell-Yan process qq ^ V (dashed line). The shaded areas are excluded by the 
experimental cuts for Q = 120 GeV. 


Since the function corrects only the <5(1 — x) piece of C^^^(a;), and it 

does not depend on the impact parameter b, its primary effect is to change 
the overall normalization of the transverse momentum distribution, but not its 
shape. In Ref. [2], a K-factor was defined as the ratio of the NLO resummed 
cross section to the LO non-resummed cross section, using the correspond¬ 
ing PDFs in the numerator and denominator. By approximating the function 
Vgg^^^ by the analogous one for Higgs production, Eq. (10), the K-factor for 
the process gg 77 was estimated to be 1.45-1.75. We can now consider 
the impact of the correct function on the K-factor. Given that the contri¬ 
bution of {Cg/s f-s){x,b; fip) constitutes less than 25% of the contribution 
of {Cg/g ® fg){x,b] lip) iu tfic ccutral rapidity region, we estimate the correct 
gg 77 X K-factor to be about 1.2-1.5. Furthermore, we can use Fig. 2 in 
Ref. [2] to find the corrected K-factor for all included subprocesses to be about 
1.3 at Q = 80 GeV and 1.6 at Q = 150 GeV. We note that the resummed 
K-factors for the gg 77 subprocess are slightly different than the fixed- 
order K-factors obtainable from Fig. 4(a) in Ref. [5]; however, this difference 
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is primarily due to the fact that the renormalization scale was chosen to be 
= Q/2 in Ref. [5] and that the different selection cuts used in that pa¬ 
per produced a kinematic enhancement of the K-factor for Q near 80 GeV. 
Of conrse, these first estimates of the corrected resnmmed K-factors can be 
further refined by repeating a detailed Monte-Carlo study as in Ref. [2]. 

Recently, it has been shown that the remaining O^ag/n'^) coefficient in the 
Sndakov factor, can also be obtained from the NLO cross section, nsing 
the universality of the real emission corrections and the general structure 
of the virtual corrections in the soft and collinear limits [12]. Following this 
argnment, we obtain 


B. 


( 2 ) 


5P. 


( 2 ) 


9/g 


+ /^o 


V, 


99^X 


+ 




12 


+ 00 111 


/iQ 

Q 


ll36 127 18 V Cl 


Pf^O^c In 


C'l, 


(In 02 


/9olnC^2, 


( 11 ) 


which is valid both for Higgs boson production and diphoton production. In 
this formnla the 5(1 — x) part of the 0{ag) splitting fnnction ioY g ^ g is 

1 ,. ^ 2 


Xp(2) _ 


(I + 3C(3)) - -N,Cp - ^N,N, , 


( 12 ) 


where ({n) is the Riemann zeta function, with ^(3) = 1.202057.... For Higgs 
boson prodnction, Eq. (11) has been corroborated by direct calcnlation from 
the NLO transverse momentnm distribntions [13]. In Fig. 2 we plot the 
coefficient functions for various processes, with the canonical choice of pa- 
rameters and Nf = 5. From this plot, we see that B^g^^ is almost exactly 
twice as large as over most of the allowed kinematic region, and both 

coefficients are considerably larger than those for the gg-initiated processes. 
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Figure 2. Comparison of the coefficient in various particle reactions. The shaded 

areas are excluded by the experimental cuts for Q = 120 GeV. 

In the standard CSS formalism, the fnnctions B and Cg/a are process-dependent, 

as seen explicitly above. Ref. [14] proposed a modihed resummation formula, 

which removes from these functions all terms associated with hard virtual 

QCD corrections to the LO process. Such hard corrections are absorbed in a 

new function H^as), so that the alternate formula for W{h,Q,v,XA,XB) is 

W(b,Q,v,XA,XB) = 

X exp |- % U(ax(p)) 1.1 ^ + B'(ax(p))l | 

{ Sclli? fiA J J 

X E [c',i. ® /a] {xa, b] f^F, 0's{j^F)) ® fb] (xb, h] /iF, asi^ip))- (13) 

a,b 

Here we can expand B\as) and Cg^^{x,b; fiF,o:s) as a series in 0:5 exactly 
as the functions B{as) and Cg/a{x^b] ^,F-,Ois), and the function 1-L{as) can 
similarly be expanded as 'H{as) = 1 + (as/Tr)” 

In this formulation, there is a “scheme-dependent” ambiguity in the dehnition 
of Cg/g, B', and H, since a change in H can be compensated by redehnitions 
of Cg^g and B'. A reasonable choice of scheme is to dehne 

W'‘Lx = h|21 + 2Aln^, (14) 

SO that C^gjg{x) vanishes for the canonical choice of parameters. In this scheme. 
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which is similar to the ‘NS resummation scheme’ of Ref. [14], we obtain 


C-^g/g{x) =S{l-x) (-Po In - Nr. In^ - P[]l(x) In 


C'l 


C'l 


a/g^ 


( 15 ) 


B. 


{ 2)1 _ 


5P. 


|^ + /3„/V,^+AAfJ (in^j -(InC^) 


( 16 ) 


as well as and ~ ^o/e(^)- The advantage of this formu- 
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lation for diphoton production is that it allows us to shift all dependence on 
the kinematical variable v from and B' into the single hard factor Ti. This 
choice makes sense physically, since this kinematical dependence is a property 
of the hard gg 77 process, rather than of soft or collinear effects. This for¬ 
mulation also makes more obvious the fact that the function affects the 

normalization, but not the shape, of the transverse momentum distribution. 
A similar modification can be made to the qq 77 X resummation formula. 


In conclusion, we have calculated the remaining unknown parts at 0{as/'n') 
in the resummed cross section for the production of photon pairs in gluon- 
gluon fusion at small Qt- We found that the approximation of the function 
Cgig{x, b; hf) in the process gg 77 by its counterpart from Higgs boson pro¬ 
duction overestimates the gg 77 X resummed K-factor by about 15 — 20%, 
and it overestimates the K-factor for the total diphoton production process 
by about 5 — 10%. We have also calculated the OicPgjTp) coefficient in 
the perturbative Sudakov factor. We predict that the impact of the coeffi¬ 
cient on the shape of transverse momentum distributions in gluon fusion 
is more substantial then in the process qq 77 , and that it will improve 
the matching of the resummed calculation with the fixed-order calculation at 
intermediate Qt ■ 
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